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Abstract 

This work presents a new approach to find the exact 
solutions for the free vibration analysis of a beam based on 
the Timoshenko type with different boundary conditions. 
The solutions are obtained by the method of Lagrange 
multipliers in which the free vibration problem is posed as a 
constrained variational problem. The Legendre orthogonal 
polynomials are used as the beam eigenfunctions. Natural 
frequencies and mode shapes of various Timoshenko beams 
are presented to demonstrate the efficiency of the 
methodology. 
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Introduction 

Beams play an important role in the creation of 
mechanical, electromechanical, and civil systems. 
Many of these systems are subjected to dynamic 
excitation. As a consequence, the exact determination 
of the natural frequencies and mode shapes of linear 
elastic beams have been studied by many researchers. 
It has been known for many years that the classical 
Euler-Bemoulli beam theory is able to predict the 
frequencies of flexural vibration of the lower modes of 
thin beams with adequate precision. The vibratory 
motion of thick beams is described by the Timoshenko 
beam theory, as they incorporate the effects of rotary 
inertia and deformation due to shear. During the past 
decades, the free vibrations of Euler-Bernoulli beams 
have received considerable attention of many 
researchers, but only few publications were devoted to 
including the effects of shear deformation and rotary 
inertia. 

The mode shape differential equation describing the 


transverse vibrations of a hanging Euler-Bernoulli 
beam under linearly varying axial force has been 
derived by Schafer (1985). Lee and Ng (1994) have 
computed the fundamental frequencies and the critical 
buckling loads of simply supported beams with 
stepped variation in thickness using two algorithms 
based on the Rayleigh-Ritz method. The first 
algorithm which has been used extensively in 
analyzing beams with non-uniform thickness, involves 
using a series of assumed functions that satisfy only 
the external boundary conditions and disregard the 
presence of the step. The second algorithm considers a 
beam with a step as two separate beams divided by 
the step. Two different sets of admissible functions 
which satisfy the respective geometric boundary 
conditions have been assumed for these two fictitious 
sub-beams. Geometric continuities at the step have 
been enforced by introducing artificial linear and 
torsional springs. 

Lee and Kes (1990) have conducted a study to 
determine the natural frequencies of non-uniform 
Euler beams resting on a non-uniform foundation with 
general elastic end restraints. The free vibration 
response of an Euler-Bemoulli beam supported by an 
intermediate elastic constraint has been studied by 
Riedel and Tan (1998) using the transfer function 
method. Rosa and Maurizi (1998) have investigated 
the influence of concentrated masses and Pasternak 
soil on free vibration of beams and gave exact 
solutions for Bernoulli-Euler beams based on the 
beam theory. A modified finite difference method has 
been presented by Chen and Zhao (2005) to simulate 
transverse vibrations of an axially moving string. 

Lin and Tsai (2007) have determined the natural 
frequencies and mode shapes of Bernoulli-Euler 
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multi-span beam carrying multiple spring-mass 
systems. An analytical solution has been presented for 
the natural frequencies, mode shapes and 
orthogonality conditions of an arbitrary system of 
Euler-Bemoulli beams interconnected by arbitrary 
joints and subject to arbitrary boundary conditions by 
Wiedemann (2007). Failla and Santini (2008) have 
addressed the eigenvalue problem of the Euler- 
Bernoulli discontinuous beams. A simulation method 
called the differential transform method (DTM) has 
been employed to predict the vibration of an Euler- 
Bernoulli beam (pipeline) resting on an elastic soil by 
Balkaya and Kaya (2009). Alim and Akkurt (2011) 
have investigated the free vibration analysis of straight 
and circular beams on elastic foundation based on the 
Timoshenko beam theory. Ordinary differential 
equations in scalar form obtained in the Laplace 
domain are solved numerically using the complementary 
functions method. 

He and Huang (1987) have used the dynamic stiffness 
method to analyze the free vibration of continuous 
Timoshenko beam. The full development and analysis 
of four models for the transversely vibrating uniform 
beam have been presented by Han et al. (1999). The 
four theories namely the Euler-Bemoulli, Rayleigh, 
shear and Timoshenko have been considered. Zhou 
(2001) has studied the free vibration of multi-span 
Timoshenko beams by the Rayleigh-Ritz method. The 
static Timoshenko beam functions have been 
developed as the trial functions in the analysis which 
are the complete solutions of transverse deflections 
and rotational angles of the beam when a series of 
static sinusoidal loads acts on the beam. 

A study of the free vibration of Timoshenko beams has 
been presented by Lee and Schultz (2004) on the basis 
of the Chebyshev pseudospectral method. Chen et al. 
(2004) have proposed a mixed method, which 
combines the state space method and the differential 
quadrature method, for bending and free vibration of 
arbitrarily thick beams resting on a Pasternak elastic 
foundation. The Laplace transform has been used to 
obtain a solution for a Timoshenko beam on an elastic 
foundation with several combinations of discrete in- 
span attachments and with several combinations of 
attachments at the boundaries by Magrab (2007). 

In the present paper, a novel approach is made to the 
problem of the free vibrations of a Timoshenko beam, 
in which the orthogonal Legendre polynomials in 
conjunction with Lagrange multipliers are used. The 
frequencies and the corresponding mode shapes for 


common types of boundary conditions are compared 
extremely well with the available solution. 

Problem Formulation 

Consider a straight Timoshenko beam of length L, a 
uniform cross-sectional area A(=b*h), the mass per unit 
length of m, the second moment of area of the cross- 
section I, Young's modulus E, and shear modulus G. It 
is assumed that the beam is made of a homogenous 
and isotropic material. 

The kinetic energy T and the strain energy U of the 
vibrating beam can be written as: 

T = — Jg |ww, 2 + mr 2 y/ 2 ^dx (1) 

U = ^\ L 0 {My/ r +Vy)dx = { Ely/ 2 + kAGy 2 ]dx (2) 

Where y represents the shear angle ( y = w i - y/ ), 
w(x,t) and y/(x,t ) are transverse displacement and the 
cross-section rotation due to the bending moment, x is 
the axial coordinate of the beam, r is the radius of 
gyration ( = *Jl / A ) and k is the beam cross sectional 
shape factor. Also M and V are the bending moment 
and shear forces, respectively. Comma denotes 
differentiation with respect to x or t. 

Applying Hamilton's principle, the governing 
equations of motion and boundary conditions are 
obtained as follows: 

mw„ - (kAG(w ■■-!//)) x = 0 ( 3 ) 

mr 2 y/ tt -kAG(w s -y/) - (Ely/ x ) x = 0 (4) 

(MSy/)' 0 = 0, (VAw) 1 , = 0 (5) 

in other words, at the ends x =0 and L, we have: 

[either M = 0 or y/ = 0 is specified, 

[either V = 0 or w=0is specified. ^ 
The equation (6) gives the boundary conditions of the 
present case. 

Analytical Solution 

In the present work, series of solutions in conjunction 
with the Lagrange multipliers are used to study the 
free vibration characteristics of the beam. The main 
advantage of the Lagrange multiplier technique is that 
the choice of the assumed displacement function is 
easy because they do not have to satisfy the boundary 
conditions of the problem. In the present study, the 
simple Legendre polynomials are chosen as 
displacement functions, and this simplifies the 
problem further since the orthogonality properties 
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lead to simple energy expression. 

Harmonic solutions for the variables w(x, t) and 
y/{x,t) are assumed as: 

h<x;/) = W(x)e wn , if/(x,t) = 'P(x)e im (7) 

in which variables W(x) and 'Fix) are the 
displacement functions and co is the circular 
frequency. As mentioned above, the displacement 
functions can be expressed in terms of the simple 
Legendre polynomials and are given by: 


i(-\) m w m =o, t(- irx=o, 

m-0 m-0 

n, n, (10a) 

I W m = °» 1^=0 

m-0 m-0 

Clamped-Hinged Beam (C-H): 

i(-ir^„=o, i(-i) mi F m =o 

m=0 m=0 

I mU j (10b) 

£w m =0, i'P m L i J (2w-4* 1 -l) = 0 

m=0 m= 1 k x =0 


W(x)= ± W m PJx), '¥(x)=t V m PJx) 


( 8 ) 


Here, Pm is the simple Legendre polynomial of degree 
m. It should be mentioned that the axial coordinate is 
transformed to the interval -1 < x < 1 by letting 
x-L/2 


We have four boundary conditions for each beam, i.e. 
two boundary conditions for two ends. These four 
boundary conditions which are not satisfied by the 
assumed series, are imposed as constraints. For four 
common boundary conditions (B.C.s), these 
constraints can be written as follows: 


Clamped-Clamped Beam (C-C): 

W{- 1) = 0, 
W( 1) = 0, 


'P(-I) = 0 
'F(l) = 0 


(9a) 


Hinged-Hinged Beam (H-H): 


m-1 

~Y 


t (-l)'X =0,1 t (-\r 2kr '{2m-4k t -1) = 0 


m= 0 


m = 1 =0 

m-1 

~Y 


(10c) 


±w m = 0, i'¥ m I (2m-4^-l) = 0 

m= 0 m=l k { =0 

Clamped- Free Beam (C-F): 


L(-iyX,=o, 

m = 0 

I m-1 

2 "< \-~2~ 


I (-1)"^ =0 


ylK I (2m-4C, -l)- I ^F,„ =0 (lOd) 

/v m=l =0 m= 0 


Z (2/77 — 4Aq -l) = 0 


m=l A:, =0 


Clamped-Hinged Beam (C-H): 


W(- 1) = 0, 
W(\) = 0, 


V P(— 1) = 0 

'P'(l) = o 


Hinged-Hinged Beam (H-H): 

JV(-1) = 0, 'P’(-l) = 0 

tR(l) = 0, 'P'(l) = 0 

Clamped- Free Beam (C-F): 

W(- 1)=0, V F(-1) = 0 

jW'( 1)-'F(1) = 0, 'F'(l) = 0 


(9b) 


(9c) 


(9d) 


in which prime denotes differentiation with respect to 
x. The boundary conditions yield linear constraints 
related to the linear combinations of the Legendre 
polynomials and subject to the degree of 
approximation, i.e. the number of polynomials 
involved. By substituting equations (8) in equation (9), 
the constraints can be rewritten as: 


Clamped-Clamped Beam (C-C): 


In above equation, we have used the following 
properties of Legendre polynomial (Gradshteyn and 
Ryzhik, 2007): 

IT] 

P;,(x)= Z (2n-4k-\)P n _ lk Jx) (n > 1) 

k=0 

in which ^ - l j signifies the integral part of (n-l)/2 
(Gradshteyn and Ryzhik, 2007). 

A variational principle is formulated based on the 
kinetic and strain energies by a procedure similar to 
one followed by Washizu (1982). This variational 
principle along with the constraint conditions is used 
to solve the vibration problem. The function to be 
extremized is given by the expression: 

4 

F = U + T - Z (Xf (constraints equation) 

where a i (/ = l-4) are the Lagrange multipliers. 
Substituting the assumed series for W(x) an dF(x) in 
equation (11) and simplifying yields: 
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m - 1 , 


n-l 

n l 

tJ 

n t 

[t 


FI 

F = T l-i Z ^ Z /I J (2»-4t 2 -\)P n _^ 4 (x)dx 

L m=l fc=0 1 ’ A 


n=l i= 


2 2 
4 io2m+l " 



ffl-1 


H 

”, 

T_ 

n t 

. 2 . 


+— tilV. Z {^k r \)P m _ lk4 (x)lW n l (2n-^-\)P n _ lii4 (x)dx 

L m=l \= 0 n=l i 2 =0 

n-1 

m = 0 «=1 £,=0 

v 2 /u/2 , 2\ v 


( 12 ) 


4 m =o 2m + 1 


l — ^ ( W m + r'P m ) - X a, (constraints equation) 

H 


The necessary extremizing conditions are given by: 

dF dF n n , 0 

= = 0 , m = 0 , 1 , 2 ,... 

SIR, d x ¥ m (13) 


Using equation (13) in conjunction with equation (12) 
results in a system of linear algebraic equations which, 
in matrix form, can be written as: 

a{w 0 ,w l ,...,w n ,'¥ 0 ,'¥ l ,...,'¥ n } T =b (14) 

in which the right hand side of equation (14) consists 
of Lagrange multipliers. Solving equation (14) for 
W m and *¥ m (m - 1,2, and substituting into the 
constraint equations (10) results in a system of 
homogenous linear algebraic equations with the 
Lagrange multipliers as unknowns. The system of 
equations is given by: 

£{a l ,a 2 ,a 3 ,a 4 } T =0 

The natural frequencies and corresponding mode 
shapes of beams can be calculated using equations (14) 
and (15). In calculating the natural frequency, the 
determinant of the coefficient matrix in equation (15) 
is computed for various values of frequency starting 
from a near zero value. Determinant change of sign 
function is identified and the corresponding value of 
frequency is the natural frequency of the beam in 
question. 


Results and Discussion 

In order to demonstrate the high accuracy of the 
present method, the convergence and comparison 
studies are carried out. Unless mentioned otherwise, 
in all of the following analysis the rectangular cross- 
sectional beams with shear correction factor k=5/6 and 


the Poisson ratio u=0.3 are considered. The first five 

I m [p 

dimensionless frequencies ( D. = J — 4 a > ) of hinged- 

hinged (H-H) and clamped-clamped (C-C) beams are 
given in Table 1. The number of terms of the Legendre 
polynomial steadily increases from 6 to 10. One can 
see that the convergence is very rapid. In general, 10 
terms of the Legendre polynomial are enough to give 
satisfactory results. 

Table 1 The convergence study on the first five dimensionless 
FREQUENCIES OF H-H AND C-C TIMOSHENKO BEAMS FOR (L/H=10) 


B.C. 

nt 

£2i 

0,2 

£^3 

04 

05 

H-H 

6 

2.1251 

5.4903 

10.1815 

21.8025 

34.7663 

7 

2.1251 

5.4467 

10.0291 

15.8141 

30.3098 

8 

2.1251 

5.4461 

9.8518 

15.4366 

22.2419 

9 

2.1251 

5.4456 

9.8425 

14.9628 

20.6510 

10 

2.1251 

5.4456 

9.8425 

14.9628 

20.6507 

C-C 

6 

0.9836 

3.7670 

8.1029 

17.3346 

37.0371 

7 

0.9836 

3.7588 

7.9795 

13.6489 

26.0859 

8 

0.9836 

3.7588 

7.9211 

13.2646 

20.2042 

9 

0.9836 

3.7588 

7.9189 

13.0420 

18.8547 

10 

0.9836 

3.7588 

7.9189 

13.0419 

18.8543 


The comparison study has been given in Table 2 for 
the first five dimensionless frequencies of Timoshenko 
beams by using the present method, the dynamic 
stiffness method (DSM) (He and Huang, 1987) and 
static Timoshenko beam function (STBF) (Lee and 
Schultz, 2004). Three types of boundary conditions: H- 
H, C-H and CC have been considered. Excellent 
agreement has been observed for all cases, which 
shows that the present method has very high accuracy. 


Table 2 The comparison study of the first five dimensionless 

FREQUENCIES OF H-H, C-H AND C-C TIMOSHENKO BEAMS FOR 
(L/H=6.667) 


B.C. 

Methods 

£2i 

£^2 

03 

04 

05 

H-H 

Present 

0.9644 

3.5195 

7.0424 

11.0702 

15.3445 

DSM 

0.9644 

3.5194 

7.0424 

11.0702 

15.3444 

STBF 

0.9644 

3.5194 

7.0424 

11.0702 

15.3444 

C-H 

Present 

1.4387 

4.1632 

7.6625 

11.5809 

15.7282 

DSM 

1.4386 

4.1632 

7.6625 

11.5807 

15.7273 

STBF 

1.4386 

4.1633 

7.6626 

11.5814 

15.7293 

C-C 

Present 

1.9814 

4.7859 

8.2462 

12.0601 

16.0888 

DSM 

1.9814 

4.7859 

8.2461 

12.0580 

16.0887 

STBF 

1.9814 

4.7860 

8.2462 

12.0604 

16.0889 


Fig. 1 shows the effect of length-to-thickness ratio 
(L/h), where L varies while h keeps constant, on the 
dimensionless fundamental frequency of the beam 
with C-C, C-H, H-H and C-F boundary conditions. It 
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is clear that for the beam with L/h<30, the variation of 
L/h has drastic effect on Q. While for higher values of 
L/h, the fundamental frequency tends to be constant, 
that is, the influence of L/h is practically negligible. As 
expected, the curve shows clearly that the smaller the 
length-to-thickness ratio is, the lower the frequency 
will be. 
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FIG. 1 EFFECT OF LENGTH-TO-THICKNESS RATIO ON THE 
DIMENSIONLESS FUNDAMENTAL NATURAL FREQUENCIES 
OF THE BEAM WITH VARIOUS BOUNDARY CONDITIONS 

Fig. 2 shows the variation of the displacement Wand 
the bending slope V F along the beam length for the 
first three modes of vibration of Timoshenko beam 
with thickness ratio L/h=10 and C-C boundary 
condition. 

Since the conventional beam theories can not involve 
the effect of Poisson's ratio, it is rather interesting to 
take a deep insight into it using the present approach. 
Table 3 gives the variation of the first three natural 
frequency parameters (Q) of H-H beams with the 
Poisson's ratio. It is shown that the natural frequency 
decreases gradually with the increasing of Poisson's 
ratio. We can see that the natural frequencies for u=0.5 
have an apparent deviation from that for u=0.1. From 
this point of view, the Poisson's ratio is of great 
significance in structural design especially for 
composite material beams. 



MODE 1 



MODE 2 



MODE 3 


FIG. 2 VARIATION OF DISPLACEMENT AND THE BENDING 
SLOPE ALONG THE BEAM LENGTH (MODES 1-3) 

W( ),T( ) 


Table 3 Effect of Poisson's ratios on the first three natural 

FREQUENCY PARAMETERS (Q) OF H-H BEAMS 


L/h 


Poisson's ratio (u) 

0.1 

0.2 

0.3 

0.4 

0.5 

10 

Qi 

0.9854 

0.9845 

0.9836 

0.9827 

0.9817 

Q 2 

2.1948 

2.1897 

2.1821 

2.1870 

2.1818 

Q 3 

3.7844 

3.7714 

3.7586 

3.7459 

3.7334 

25 

Qi 

0.9978 

0.9977 

0.9975 

0.9974 

0.9972 

Q2 

3.9625 

3.9600 

3.9576 

3.9551 

3.9526 

Q3 

8.8136 

8.8017 

8.7900 

8.7782 

8.7665 


Conclusions 

The free vibration of the Timoshenko beams is 
investigated using an assumed series solution in 
conjunction with Lagrange multipliers. It is observed 
that the present method is a computationally efficient 
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tool in predicting the natural frequencies of the beams. 
This method is particularly attractive because of the 
ease with which one can choose the generalized 
displacement functions. This fact is shown by 
choosing Legendre Polynomials whose orthogonal 
properties simplify energy expression considerably. 
The natural frequencies of the Timoshenko beam 
obtained by this method compare extremely well with 
the available exact solution. It should be mentioned 
that by applying this method, the convergence is very 
rapid. 
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